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I. 



INTRODUCTION 



An important eleaent in Naval Warfare Defence is 
the development of a surveillance system capable of 
monitoring the appearance and progress of any potential 
threat. The infrared (IR) spectral region shows good 
promise for passive detection systems able to detect a 
variety of targets at different altitudes. Such a system 
is the Infrared Search and Target Designation (IRSTD) 
System (AN/ SAR-8) currently in engineering development. 

In support of the development of such a system, it 
is necessary to fully specify and understand the IR 
radiance from the earth - atmosphere scene within a 
given f ield-of -view (FOV). Detecting a target by its 
radiation requires that not only the target intensity but 
also the spectral and spatial structure of the background 
be fully understood. 

Important inputs to the design of successful 
detection systems are recorded magnitudes and variations 
in the radiance from various scenes. This information 
will help in developing data processing algorithms to 
distinguish radiation signatures of clouds and other 
surfaces from those of the target. 

In the context of this THESIS , data on cloud and 
background radiance were collected and analyzed for the 
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purpose of later development of algorithms that will aid 
the AN/SAR-8 to ignore false targets generated by 
background , atmospheric or non atmospheric, IR radiation. 

These data were collected using the M AGA 780 
Thermovision" , a liquid nitrogen cooled thermographic device 
able to detect IR radiation in the Middle IR, called 
hereafter Short Wave - SW (3-5.6 micrometers) , and in the 
Far IR , called hereafter Long Wave - LW (8-14 micrometers). 

The present work consists of eight chapters following 
the introduction . The first chapter describes the 
problem and the expected output of the work. In 
the second chapter a theoretical overview is 
presented , along with data from the literature , of 
clouds and different backgrounds. The third chapter 
describes how the LOWTRAN propagation / radiance code 
was used to give transmittance of the atmosphere under 
the various circumstances of the measurements. The 
fourth chapter lists the capabilities of the AGA 
equipment and the mode of operation used . The fifth 
chapter describes how the AGA was calibrated 
empirically . In the sixth chapter the method which was 
actually used is presented both mathematically and 
analytically . In the seventh chapter measured cloud and 
background signatures are listed, analyzed and presented 
as final results. The last chapter consists of the 
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conclusioDs that the results lead 



to and 



suggestions 



for possible inproveaents and future work in the area. 



The results 


obtained 


are not 


expected 


to 


have 


significant quantitative 


iaportance 


, due 


to 


the 


wide variety 


of possible background radiators 


and 


aeteoro logical 


conditions , 


but it is 


hoped 


that 


they 


will provide a 


starting 


point for 


further 


work 


on 


the subject. 
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II. BACKGROUND TO PROBLEM 



A. GBNERAL 

All objects in our environment , either cold or hot, 
radiate IR energy detectable by the IRSTD system. 

The IRSTD , as mentioned before , should be capable of 
determining whether an observed IR source is (1) an airplane/ 
missile , or (2) an unwanted source such as clouds, birds, 
or shore features. 

There are several atmospheric and terrain conditions 
affecting the translation of emitted radiation into 
meaningful signal. Two of the most important factors are the 
atmospheric background and the atmospheric absorption. 

Atmospheric background consists mainly of clouds, rain, 
snow, clear sky etc. Terrestrial background consists of 
land, mountains, hills, cities, buildings etc. 

Here we will deal with cloud and building 
signatures. The radiance of a selection of such sources was 
measured using the M AGA 780” system. 

The equipment used is designed to provide the 

temperature of an object if its emissivity is known 
or vice versa. An empirical calibration of the ”AGA 

780” along with the application of some simple 
mathematical formulae was needed in order to obtain 
radiance measurements • 
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The second factor mentioned above , the atmospheric 
absorption, has been dealt with by using the LOWTRAN code for 
propagation of IR radiation. Specifically the LOWTRAN 6 
computer program was used to provide attenuation factors 
under different atmospheric conditions. 

B. EFFECTS OF BACKGROUND CLUTTER IN IMAGING 

Background clutter is an important factor 
limiting the performance of thermal imaging equipment . 
Clutter emission from clouds , birds and terrestrial 
objects generates a sensor response which may either 
mask the signal from a target or appear as a 
false target . 

It is impossible to predict how large spatially 



or 


how intense a 


clutter 


signal 


from 


background 


would 


be. 


Background clutter may 


consist 


of 


either 


emission 


or 


scattering from 


smal 1 


objects 


(eg. 


birds ) 


as 


well 


as 


from large 


objects 


(eg. clouds) 


[Ref. 1] 


• 






Targets and 


clutter 


will 


be 


sensed 


by 


an 



imaging system as having a small temperature 

difference from the background, AT. This temperature 

difference can be expressed as radiation contrast C , 

K 

as described by Lloyd [Ref. 2]. 



R 



. w i_:_ w b. 

W T + W B 



(II-l) 



W T = target radiant emittance ( W/m ) 
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2 

W n = background radiant enittance ( W/a ) 

D 

If the total radiant eiittance W(T) is 

W(T) = or T 4 (Stefan’s law) (H-2) , then 

*W(T) „ 

= 4 or T J ( II-3) 

91 

or, AW(T) = 4 or T 4 AT (II-4) 



Now 



For 



R 



W 



B 



W T ♦ w B 

W(T) = or T 4 



(II-l) 
( II-5) 



and W T 



C 



R 



C 



R 



3W(T) 

W(T+ AT) = W(T) + AT (II-6) 

5T 



AW(T) 



2W( T) + AW ( T) 

4 or T 3 AT 

2 or T 4 + 4 or T 3 AT 
2 AT 

or, 

T + 2 AT 

2 AT 

, for small AT. 

T 



(H-7) 



( 1 1-8 ) 



For typical scenes AT is of order 1°K so this 
approxiaation nay be justified. 



C. BACKGROUND SUPPRESSION SCHEMES 

A filtering systen ( subtraction scheme ) will be 
needed to suppress the background signal. This could 
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be an automatic background suppression system that 
will sense the average background signal and subtract 
it . 

A problem could arise here when the temperature 
of the target and its background are nearly the 
same ( zero contrast ) . This will happen twice each 
day . 



Hudson [Hef. 3] gives a very 


good exaipie 


showing that last point. A truck is 


parked in an 


open field so that an infrared system 


can view it . 


Observations are made over a period of 


24 hours and 


variations are found in the contrast 


between the 


vehicle and its background . In the 


afternoon the 


truck has been heated enough by the 


sun to be 



warmer than the background , giving a positive contrast. 

During the early evening hours the vehicle , because 

of its large thermal capacity , cools more slowly than 
the background and the contrast is even greater ( more 
positive ) . At night the truck cools more rapidly , 
passes the point of equal temperature to the 
background and the contrast passes through zero to 



become negative. As the sun rises again 


the background 


warms more rapidly than the truck and 


the contrast 


becomes more negative . Later in the 


morning the 


heating of the sun will be sufficient 


to cause a 



period of zero contrast and then a positive one. That 
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shows that ” for many targets and background 
combinations , there are two intervals in any 24-hour 
period during which the target cannot be detected 
because there is insufficient radiation contrast 
between it and its background " . 

There are no good methods of eliminating this 
effect ; the best thing to do would be to wait until 
the period of "washout" passes. 

D. CLUTTER SUPPRESSION SCHEMES 

Assuming that targets and clutter will have 
different temperature from the background and that the 
background suppression is successful , the problem of 
eliminating clutter still remains. 



Several unwanted 


signals 


that will 


have 


small AT 


above 


the background 


will be 


present . 


These 


signals 


can be 


from birds 


, clouds 


, buildings 


etc 


. . They 


will all appear as 


threshhold 


crossings 


exactly as a 


target 


could appear. 










To 


avoid confusion caused by 


these 


signals 


several 


nethods can 


be used . 


The use 


of 


each one 



will depend on the task that the system will 
perform ( early warning, detection , tracking etc. ) and 
on the kind of clutter expected ( small or large , 
intense or not ) . 
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One method could be the use of a rotating 
reticle , in the equipment’s image plane , to 
achieve space filtering . This method is described by 
Hudson [Ref. 3] and can be used for discrimination 

of small targets against large clutter like clouds in 
seeker or tracker systems. 

Other methods will be dependent mainly on the 

size of the target to be detected . 

For example for a small target we could divide 



the 


scan 


electronically 


into 


small windows 


( each 


one 


being 


equal 


in 


dimension 


to 


one Instantaneous 


field 


~of~ 


view 


) and 


test 


if the 


target appears in 


more 


than 


two 


neighbouring 


windows . 


In 


figure 2.1 signals A 


and 



B will be classified as targets and signal C will be 
classified as clutter and discarded. This scheme will 
not be effective against small false targets such 
as birds . 

For large targets (eg. ships) a useful logic 

method could be to check if the target appears in 
three out of four successive windows. In figure 2.2 
signals A and B will be classified as targets but 

signal C will be classified as clutter. This system 
may perhaps not detect a small target (however , a 
missile may be detected since its booster is heating 

the air surrounding the exhaust gas , and will give 
a large target ) . 
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that will detect small targets 
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Fig. 2.2 Clutter suppression scheme for system 

that will detect large targets. 



Other possible Methods could be software 
of the noveient of the targets (eliminating for 
all signals not showing any Movement at the 
discarding ships that will be stationary) , or 

of signature patterns and statistical comparison 
known signatures ( although this might be 

difficult due to the wide variety of 

signatures ) . 



analysis 
example 
risk of 
analysis 
with 
very 
clutter 



E. CONCLUSION 

What is obvious from the above is that it will 

be useful to know what amount of thermal radiation 

is emitted by different kinds of clutter like clouds, 
buildings etc, . The purpose of this THESIS is to 
provide and analyze thermal radiation values from such 
clutter. 
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Ill* CLOUDS * BACKGROUNDS AND THEIR SIGNIATURBS 
A. CLOUDS 

Clouds are one of the sain causes of undesirable IH 
radiation . Cloud cover in lost of the regions of the world 
is a condition that exists about one half of the tine. Figure 
3.1 shows the frequency of overcast skies (curve 1) and 
the total frequency of overcast and partly cloudy skies 
(curve 2) along the 20° W Meridian, as a function of 
Latitude [Hef 5]. 



O 




LATITUDE 



Fig. 3.1 : Latitude and dependence of frequencies of 

various kinds of weather along the 20° W 
■eridian.(l) frequency of overcast skies 
(8-10 on a 10 fold scale for cloudiness) 

(2) frequency of overcast and partly 
cloudy (3-7) skies . Region above curve 2 
characterizes frequency of clear skies. [Ref . 5] 



22 



Clouds may be classified into three Major categories 
according to the height at which they occur. 

(i) High level clouds. 

(ii) Middle level clouds. 

(iii) Low-level clouds. 

The following table shows those categories* the type of 
clouds in each and typical heights at which they occur. 



Cloud 

level 


Type of 
Clouds 


Altitude at which they 
occur at Midlatitude. 
(Avg. all seasons) 




Cirrus 






High 


Cirrocumulus 

Cirrostratus 


5-13 


km 


Middle 


Altocumulus 


2-7 


km 


Low 


Stratocumulus 


Earth* 


s surface 




Stratus 


to 2 


km 



Figures 3.2 and 3.3 both from [Ref. 61 show 
spectral radiance from two types of clouds , Cumulus and 
Cirrus respectively. 

Discussion of the particular behavior of each category 
is beyond the task of this THESIS but in general we can state 
that the clouds perform three things of interest in the 
present context: 

(i) They reflect radiation originating from 
terrestrial objects. 
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Figure 3.2 : The spectral radiance of the underside 

of a dark cuaulua cloud . [Ref. 6J 




Figure 3.3 : The spectral radiance of sky covered 
with cirrus clouds at several angles 
of elevation. [Ref. 6] 
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(ii) They scatter solar radiation strongly , and 

(iii) They emit IH radiation due to their 
own temperature. 

Figure 3.4 [Ref. 10] shows how solar thermal 
radiation components reach the surface of the earth and 
then how they are reemitted back. There we see that of 
the incoming thermal radiation , 3% is absorbed by 

clouds and 20* is reflected by clouds back to space. 
On the other hand emission by clouds is responsible 
for 26X of the outgoing radiation to space. 

1 . Reflection of IR Radiation by Clouds 

Clouds reflect IR radiation from terrestrial objects. 
Typically the reflectance of a middle layer cloud ranges from 
0.75 to 0.1 in the spectral range 0.2 to 4 pm (fig. 
3.5 ) [Ref. 6]. 

The reflectance of a cloud increases with the cloud 
thickness, and with decreasing mean free path of a Light ray. 
It also varies with the solar zenith angle (fig. 3.6) 

[Ref. 6] . 

2 . Scattering of IR Radiation by Clouds 

The scattering process is the process by which 
particles absorb energy from an incident electromagnetic 

wave , from some direction , and reradiate that energy into 

the total solid angle centered at the particle. Scat- 

tering is a function of particle size, density and 
refractive index and also wavelength of incident light. 
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Figure 3.4 : A Buiiary of the solar and thermal 

radiation conponents of the Earth’s 
Radiation budget. Numbers are percentages. 
[Ref. 10 ] 




Clouds. Data are directional reflectance of a middle laver cloud. 

• •••••••• Winter Snow and Ice. Data are directional reflectance of dry snow. 

............ Summer Ice. Data are directional reflectance of summer Arctic ice. 

• - • - • - Soil and Rocks. Data represent the average value of the bidirectional 

reflectance, o x (45°. 0. 0. 0), of gravel, wet clay, dry clay, 
tuff bedrock, and sandy loam. 

— . — _ Vegetation. Data represent the average value of directional reflectance 
ot many types ol vegetation (from the EK1M Data File). 



Figure 3.5 



Spectral diffuse reflectance of 
ataosphere constituents. [Ref. 6] 



Earth- 
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Cloud Reflectance 



Figure 3.6 : The variation in reflectance as a 

function of cloud thickness parameter 
hL for different solar zenith angles. 
The cloud thickness is h and the 

■ean free path of a light ray is 
L . [Ref. 6] 
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Water clouds contain water-droplets and ice-crystals 

with radii from 1 fin to about 10 fin and concentrations of 
3 

300-10 per cn . (for comparison purposes the corresponding 

2 4 —2 -3 3 

values for raindrops are 10 —10 pm and 10 - 10 per cm ) . 

Figure 3.7 [Ref. 11] shows the water droplet size 

distribution as a function of drop radius for various cloud 

cases. The highest clouds are composed of large ice crystals 

of radii greater than 50 pm. The lower clouds can be assumed 

to be at above freezing temperatures and therefore 

water droplet clouds. Figure 3.8 [Ref. 10] shows the 

composition of clouds vs. temperature. 

Under those considerations Mie scattering (particle 

size comparable to wavelength) is dominant assuming that 

incident solar radiation at the clouds has a peak at about 

0.5-1 pm and land radiation has a peak at about 10 pm. 

Rayleigh scattering (particle size much less than wavelength) 

is of less importance in connection with water clouds. 

We can say that the solar scattering through 

clouds occurs at wavelengths shorter than 3 pm (Near 

IR) . Middle IR is less affected by particle scattering 

and the Far IR is barely affected at all. 

Near IR radiation exhibits strong forward 

scattering in clouds, but for a heavy overcast sky 

multiple scattering reduces this effect. 

Finally , it is clear that the relative positions of 

sun, observer and cloud cover are of importance in 
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Figure 3.7 : Drop size distributions n(r) for various 
cloud cases. [Ref. 11] 

A) stratus base 

B) stratus top 

C) stratocumulus base 

D) stratocumulus top 
R) ninbocumulus base 
F) niabocuaulus top 
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0 



f. 




Figure 3.8 



Phase composition of clouds 
temperature . [Ref. 10] 




Figure 3.9 



Note 



Vertical temperature distribution for 

various models. [Ref. 10] 

1 and 2 Cloudless atmosphere night and 

day 

3 and 4 Cloud, night without and with 

scattering 

5 Cloud day 

6 Cloud night (selective absorber) 
In curves 3 and 4 the cloud is assumed 
to be graybody. 



J 
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determining the amount of IR radiation (due to 

scattering) received through clouds. 

3 . Emission of IR Radiation by Clouds 

Like all other bodies, clouds emit IR radiation due 
to their temperatures • Usually the temperature of a 
cloud differs from the surrounding air by not more than 

some tenths of a degree. This difference reaches a 
maximum of one degree for dense cumulus clouds [Ref. 

10]. Under those considerations the temperature of a 

low cloud is about 10° C below ambient temperature at 
the sea surface. Also it is obvious that a cloud’s 

temperature decreases with height. Figure 3.9 [Ref. 10] 
shows the vertical temperature distribution for various 
models of clouds . 

Thick clouds ( thickness>100m) ie. many free paths 
in thickness and clouds with content of water droplets (or 

3 

ice crystals) greater than 0.5 g/cm are good blackbodies 
[Ref. 5]. Since their temperature is of the order of 300°K 
they emit mainly in the 8-13 pm band. Figure 3.10 [Ref. 
5] shows the spectral emissivity of cloud layers of 
thickness Az and figure 3.11 [Ref. 5] shows how the 

integral emissivity , transmittance and absorptance varies 
with Az . It is clear from that last figure that for 

Az > 60m r > 0 f and 

Az > 100m € > 1 



32 




Figure 3.10 



Spectral emissivity of cloud layers 
thickness Az . [Ref. 5] 




Figure 3.11 



Integral emittance (E), transmittance (T) 
and absorptance (A) curves. [Ref. 5] 
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Figure 3.12 [Ref. 6] shows the zenith sky 
spectral radiance from a cloud. In this figure we note 
two things: 

(a) The large variation of the cloud spectral 
radiance for different anbient air temperature. 

(b) That neglecting the absorption band centered 
at about 10 pm the cloud spectrum approaches very 
closely the blackbody curves of the ambient air 
temperature thus making the latter a very good estimate 
of the cloud’s own temperature. 

B. BACKGROUNDS 

For our purpose we may consider two general categories 
of background (surface backgrounds). 

1 . Marine Backgrounds 

The surface of the sea emits its own thermal radiation 
and also reflects incident IR radiation. The major factors 
that determine the amount of IR radiation due to a marine 
background are (1) the surface temperature of the water (2) 
the optical properties of the water (3) the amount of 
agitation of the water (4) the material properties of the 
bottom of the sea [Ref. 6] . 

For our purpose, it has been assumed that marine 
backgrounds won’t produce any effect by reflecting IR 
radiation since the reflectance of the water is less than 
0.05 for angles of incidence of 50° or less as can be 
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Figure 




. 12 



Zenith sky spectral 
large variation 
temperature. [Ref. 6] 



radiance showing 
with ambient 



the 

air 
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seen in figure 3.13 [Hef. 6]. Rough sea May change that. 
For a Beaufort 4 wind roughened sea , for the sun at 
the horizon , the reflectance reaches a maximum of about 
0.20 . But for the sun at an elevation of 30° or more 
the reflectance will be less than 0.05. (Figure 3.14 
[Ref. 6]) 

What remains is the emission of IR radiation by the 
sea. The emissivity of an average rough sea is about 0.8 but 
the amount of IR radiation actually emitted depends on the 
water temperature. This temperature is near 0°C (about 4°C) 
in the arctic regions and 29°C near the equator. 

Currents may produce anomalies of several degrees Celcius 
as warm water flows into cold water or vice versa. 

2 . Terrestrial Backgrounds 

Terrestrial objects have a wide variety of emis- 
sivities and temperatures. Some of them also are not opaque. 
Obviously the radiance of a sun heated object falls 

during the night when heating from the sun is absent. 
Figure 3.15 [Ref. 6] shows the variation in the 10 pm 
radiance of selected backgrounds during a 24 hour time 
period. 

Most terrestrial objects are not blackbodies. This 
can be seen in figure 3.16 taken from [Ref. 6], which 

represents the radiance of a city on a plain (Colorado 

Springs) as viewed from a mountain (Pikes Peak) at a 
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Figure 3 . 13 



Reflectance and eaissivity of water ( 2 

to 15 pa average ) versus angle of 

incidence. (Ref. 6] 




Figure 3 . 14 



Reflection of 
flat surface 



solar 



radiation 
and a 



f roa a 

surface 



( * s =o ) 

roughened by a beaufort 4 wind (or =0.2) 
The lower and upper branches or the 
curve aarked 9 =0.2 represent two 
assuaptions regarding the effect of 
aultiple reflection. True values are 
expected to lie between the indicated 
liaits. [Ref. 6] 
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Figure 3.15 : Diurnal variation in the 10 fin Radiance 
of selected backgrounds . j 

SS : sunset; BNT : end of nautical twilight; j 
SR : sunrise; BNT : beginning of nautical I 
twilight. [Ref. 6] I 




Figure 3.16 : Radiance of an urban area and clear 

zenith sky (Colorado Springs fron Pikes 
Peak). [Ref. 6] 
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IV. LOWTRAW PROPAGATION / RADIANCE CODE 



A. MODELS FOR PROPAGATION OF IR RADIATION THROUGH 

THE ATMOSPHERE 

For aany applications , and specifically for the 
performance of military systems in imaging , a model 

which will provide accurate transmission and modification 
of IR radiation through an atmospheric path , and also 
the fluctuation in these effects is needed. 

Two such models are available , the High resolution 

model developed by the AFCRL and the Low resolution 

model called LOWTRAN. 

The HITRAN High Resolution models such as LASER 

and FASCODE are direct computation models which provide 
very high resolution (ie. single frequency) results. This 
is necessary to predict atmospheric behavior of laser 

radiation , but not for thermal ( broad band ) sources 

such as considered here. 

B. LOWTRAN CODE 

LOWTRAN has been developed as a curve-fitting program 
matching empirical or precomputed attenuation data . This 
empirical data is averaged over some specific spectral 
interval (20 cm ^). As a consequence of the averaging 

over finite spectral interval, many of the rapid 
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spectral fluctuations of atmospheric transmittance will 
vanish. 

Thus the LOWTRAN code can be considered as a 

fast but not very accurate computer program. 

L0WTRAN6 , the latest version of LOWTRAN, is a program 
designed to calculate transmittance and/or radiance of 
the atmosphere for a specific set of geometrical and 

weather conditions , in the spectral range 0.25 to 28.57 
micrometers ( 350 to 40,000 cm at 20 cm ^ intervals (ie. 

20 cm ^ spectral resolution). 



C. INPUTS TO LOWTRAN 6 CODB 

To run LOWTRAN 6 one must specify four sets of 

input data ( cards ) . 

In Cards 1 & 2 one of the six 

atmospheric models or a user defined 
meteorological conditions ( model atmosphere ) 
specified. 

In Card 3 are specified the 

geometrical path parameters of the problem. 

In Card 4 finally are specified the 

range over which data are required and the frequency 
increment at which the data are to be printed out . 



avai lable 
set of 
must be 

different 

spectral 
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D. INPUT DATA USED 

This model has been used in this study to predict 
sky radiance and correct other radiance values for 

atmospheric effects. 

For Card 1 the Midlatitude Winter model was used. 

For Card 2 the Maritime extinction ( 23km 
visibility ) model was used most of the time. 

For Card 4 the wavelength bands were the 3 to 

5.6 pm ( 3340 to 1780 cm * ) and the 8 to 14 pm ( 1250 
to 710 cm ^ ) windows with 20 cm ^ spectral intervals. 

Otherwise the rest of the data was dependent on 

the specific conditions of each measurement. 



E. RESULTS OBTAINED 

The program was used in its radiance mode . The 

average transmittance of the atmosphere and the 
radiance of the clear sky were recorded. 

These results appear in the Table 5 of chapter 
eight along with the radiance values provided by the 
program . 
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V. EQUIPMENT 



A. THE AGA THERMOVISION SYSTEM 

"AGA Thermovision" was first introduced in 1965. The 
780 System is a fourth generation development of a system 
which combines real time IR scanning with thermal 
measurement capability. By thermal measurement is meant 
measurement of the temperature of an object. Chapter 7 
will analyze how the AGA was used to measure 
radiance. 

The ” AGA Thermovision " 780 owned by the Naval 
Postgraduate School (Electro-Optics Laboratory , Department 
of Physics) consists of a dual scanner : one 
channel for the 3 - 5.6 pm band (SW) and one for the 
8-14 pm band (LW). It uses a two dimensional object- 
plane scan. Figure 5.1 shows the different components 
of the AGA system. 

B. GENERAL DESCRIPTION OF THE "AGA THERMOVISION" 780 

1 . Optics and Scan Mechanism 

In the AGA Thermovision 780 , the scan is performed by 
two refractive prisms with eight sides. Horizontal scan 
lines are interlaced (with 4:1 interlace) to provide a 
complete frame of 400 lines. A mask at the display of 
the system allows 280 lines of the 400 to be viewed. 



43 



AGA Thermovision® 




LW BLACK/WHITE MONITOR CHASSIS 



Figure 



5.1 



Components 

[Ref. 7] 



of the 



AGA Thermovision 



780. 
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The lenses of the SW system are constructed from 
silicon, and those for the LW system are constructed 
from germanium* 

From the front of the scan head the aperture stop and 
the bandpass filter can be selected. For our purpose the 
non-filter position and the f/1.8 aperture stop have 
been used. 

2 . Detectors and Cooling 

There are two detectors , both cooled to 77°K using 
liquid nitrogen. For SW applications an Indium-Ant imonide 
(InSb) and for LW applications a Mercury-Cadmium-Tel lur ide 
(HgCdTe) detector is used. 

3 . Display 

A black and white display for each band gives 
an image of the object viewed. The brighter the 
object appears on the display, the more IR radiation 
it emits. 

4 . Technical Data 

Technical data for the "AGA 780" appear in Table 1. 

C . MEASUREMENT TECHNIQUES 

The " AGA Thermovision 780" measures , as mentioned 
before, IR radiation within two bands (windows) or 
spectral ranges, the 3-5.6 (im window and the 8-14 pm 
window. The radiation received by the equipment includes 
radiation reflected by terrestrial objects in the line of 
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TABLE 1 



AGA SENSOR TECHNICAL INFORMATION fRef. 131 



Performance characteristics 



Spectral Region_£pm) 
Frame rate (sec_^) 
Field rate (sec ) 
Interlace 

Replaceable Fore-Optics 

FOV (deg) , AZ 

EL 

IFOV (mrad) , AZ 

EL 

NEAT at 22°C ( °C) 

MDT ( °C) 

Dynamic Range ( °C) 
Optical Data 



Detector and Cooler Data 



3 to 5.6 and/or 8 to 14 (1) 

6.25 

25 

4 : 1 

F0V2 (2) 

7 

7 

1.1 

1.1 

0.12 

< 0.1 

-20 to 900 

F0V2 (2) 

24 
5.5 
9.9 
1.8 



2 

Effective Aperture Area (cm ) 
Diameter of Aperture (cm) 
Effective Focal Length (cm) 
f/ number 



Type and Material 
Number of Elements 
Peak Wavelength (pm) 



InSb or HgCdTe 
1 

5 or 10 



Notes: (1) According to [Ref. 13] the SW window is 2 to 5.6 

pm but [Ref. 7] gives 3 to 5.6 pm. 

(2) There are four other possible FOV’s (3. 5X3. 5, 
12X12,20X20,40X40) but data is presented only for 
the one used. 
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sight and also is affected by atmospheric absorption , 
thus making the relation to the temperature of the objects 
observed a non linear one* 

The equipment measures the received IR radiation 
in " Isotherm Units ( IU ) ## , an arbitrary unit of 
measurement* The value obtained has a linear relationship 
with received photon flux (photons/second). 

If object temperature is desired , the non-linear 
relation of IU and Temperature ( C) is given by 
calibration curves (SW figures 5.2 to 5.5 , LW figures 
5.6 to 5.9 ). In these curves the object surface 
emissivity has been assumed equal to one ( € =1 , 
object: black body ), and no atmospheric absorption has 
been considered. If the object is not a black body , or 
atmospheric absorption has to be taken into account, 
corrections have to be applied ( See section V.E : 
Corrections of measurements for non-ideal situation). 

There are two basic methods of thermal measurement 
using the AGA 780 : 

(i) The Direct measurement where the equipment is 
used to obtain directly the temperature of the object 
without any reference to any other source of radiation. 

(ii) The Relative measurement where the equipment 
is used to obtain temperature by comparing the radiation 
received from the observed object and the radiation 
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Figure 




AGA Infrared Systems AB 



INDIVIDUAL CALIBRATION OF 780 DUAL BBAR 



SERIAL NUMBERS 



SCANNER 

DETECTOR 

FILTER 

LENS 



SUDB 4011 
G 1739 
NOF 

7 3009 



CALIBRATION CONDITIONS: 



AMBIENT TEMPERATURE 
RELATIVE HUMIDITY 
OBJECTIVE DISTANCE 
CALIBRATION DATE 
OPERATOR 



23 C 
SO y. 

1.0 M 

82-09-17 
B. A 



CALIBRATION CURVE CONSTANTS 



APERTURE 


A 


B 


C 


i .9 


911203 


3117.74 


1 .00 


2.9 


449342 


3127. SO 


1. 00 


3.6 


267514 


3161 . 13 


1.00 


5. i 


1619S8 


3182.86 


1 . 00 


7.2 


96669 


3220 .95 


1. 00 


10 


65232 


3242 . 19 


1 .00 


14 


38173 


3253.91 


1 . 00 


20 


30288 


3405.76 


1 . 00 



.2 : SW calibration curves for the AGA 780. [Ref. 7] 

The curves that follow are of the fori 

A 




A, B, C are the empirical constants above. 
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uauorauon curves 






AGA Infrared Systems AB 



ISOTHCflM UNllc 8 




Scanner • S¥DB -4311 Date : A3 17 Operator? B*A 

Detector - G 1739 Ambient temp ; 23 £ 

Filter: Rel humidity: % 



I 



Figure 5.3 : SW calibration curves for the AGA 780. 

Temperature range Fro* — 20°C to 50°C. 
[Ref. 7] 
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Scanner i SWB ^11 Data: 32~S- 17 Operator * tLA 

Datacto n S 1738 Ambient temp* 23 ]5 

Fitter * WPP Ret humidity * ^ ^ 



Figure 5.4 : SW calibration curves For the AGA 780. 

Teaperature range froa -20°C to 150°C. 
[Ref. 7] 
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AGA Infrared Systems AB 



INDIVIDUAL CALIBRATION OF 780 DUAL BBAR 



SERIAL NUMBERS 



SCANNER 

DETECTOR 

FILTER 

LENS 



LUDB 4011 
G 1739 
NOF 

7 3105 



CALIBRATION CONDITIONS: 



AMBIENT TEMPERATURE : 24 C 

RELATIVE HUMIDITY : 55 X 

OBJECTIVE DISTANCE : 1.0 n 
CALIBRATION DATE : 82-09-27 

OPERATOR : B.A 



CALIBRATION CURVE CONSTANTS 



Figure 5*6 





APERTURE 


A 


B 


c 




1.8 


-3581 


i5Q6 . 49 


-0 .436 




2.5 


-4060 


1569.30 


-0.759 




3.6 


-6514 


1629 .70 


-1.824 




5.1 


5420 


1610 .70 


2.796 




7.2 


1123 


1606,54 


1 . 098 




10 


584 


1604.29 


0.885 




1 4 


306 


16i 0 .60 


0.773 




20 


195 


1650 .51 


0.767 


LW 


calibration curves for the 


AGA 780. [Ref. 7] 


The curves 


that 


follow are 


of the for* 






A 








I = 

C 


B/T 

e - 


1 




A, 


B, C are 


the 


empirical 


constants above. 
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Calibration curves 




AGA Infrared Systems A8 

I* 3 soTwew* uwrrs 




Scanner * UflB 4flll Daf t B2~39“27 Operator * 8- A 

Oetector t G 1739 Amb w nt temp t 24 £ 

Piter: ^ Ret humidity* ® % 



Figure 5.7 : LW calibration cuvea for the AGA 780. 

Temperature range from -20°C to 50°C. 
[Ref. 7] 
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Calibration curves 




AGA Infrared Systems AB 



H OTH CT W UNITS 




Scanner : UflB jgll Date : 32-t5-27 Operator : 3. A 

Detector : S 1739 AmOterrt tempt 24 £ 

Piter: Rei hurmdtty » ^ j» 



Figure 5.8 : LW calibration curves for the AGA 780. 

Temperature range from -20°C to 150°C. 
[Hef. 7] 
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AGA Infrared Systems AB 
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Figure 5.9 : LW cal ibrat ion 




received from an external reference source of known 

temperature and emissivity. 

For our purpose the Direct measurement method was 
used mainly because of the lack of an "ideal reference", 
ie. an object of known temperature, having the same 

emissivity and situated near the observed object. For 
example it was (almost) impossible to have an "ideal 

reference" near a cloud. 

D. DIRECT MEASUREMENT PROCEDURE 

Assuming that the emissivity of the object and 

the transmittance of the atmosphere are equal to one, 

€ =1 

o , s 



T = 1 

a 



the 


following procedure 


was 


used ( figure 5.10 and 5. 


ii 


show 


the 


equipment and 


the 


procedure 


respectively ): 






(i) 


The "THERMAL 


LEVEL 


control 


knob was set 


to 



some specific value "L" , to obtain a satisfactory picture. 

(ii) The "THERMAL RANGE" control knob was set to 
improve the picture obtained above. 

( iii ) The "ISOTHERM LEVEL 1" control knob was set at 
some value "i" in order to brighten up the object of 
interest . 

(iiii) The values "L" and "i" were added to obtain 
the radiance I in Isotherm Units. 
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12. 
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7. 
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7a. 
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14. 
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15. 


ISOTHERM WIDTH 






16. 


ISOTHERM SCALE 



Figure 5.10 



Front face of the 
■onitor showing the 
[Ref. 7] 



AG A black 
different 



and white 
controls . 
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Figure 5.11 : Direct leasureient technique. [Ref. 7] 





E. CORRECTION OF MEASUREMENTS FOR NON-IDEAL SITUATION 

Many factors affect the radiation received by the AGA 



780 froa an 


object . 


The most 


important 


ones 


and 


the 


actions taken. 


either 


to correct 


then or 


just 


to 


take 



them into account , are listed below. 

(i) Emissivi ty . 

Most of the objects have an emissivity less than 
one ( € <1) , often non-uniform over the spectral range 
( € = €(X)<1). Typically the emissivity of a cloud ranges 
from 0.2 to 0.8 but for clouds of thickness greater 
than 100 m the emissivity can be assumed to be close 
to unity, (refer to III-B-3) 

(ii) Object reflections 

An object will emit not only its own radiation 

but also it will reflect radiation from other objects 
(surroundings). If the object to be measured does not 
transmit IR radiation its reflectivity will be 
P = 1 - € ( V— 2 ) 

So , according to the assumption made in the above a 
cloud’s reflectivity will be less than 0.1 . 

(iii) Object size 

If the object does not cover the solid angle 
subtense of the detector then radiation will be 
received from the background. In most of the measurements 
made here an effort was made to have the object more 
or less covering the entire f ield-of-view . 
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(iiii) Atmospheric absorption 

Water vapor or ice (HgO), Carbon dioxide (COg) and also 

some other constituents of the atmosphere absorb the 

IR radiation. Thus IR radiation received by the system 

will be less than that emitted by the object observed. 

The ratio of these two is the transmittance of the 

atmosphere , r , and is provided by the LOWTRAN code. 
& 



60 



VI. EMPIRICAL CALIBRATION OF THE AGA 
A. GENERAL 



An empirical calibration of the "AGA Thermovision 

780" was made in order to determine if the calibration 
curves provided by the manufacturer were accurate 
enough under the circumstances in which the equipment 
was used. 

The calibration of the AGA consisted of a two 
scale experiment : 

(i) A near field experiment and, 

(ii) A far field experiment. 

B. NEAR FIELD EXPERIMENT 



The thermal radiance of an aluminium can filled 
with water was measured. The can was painted black 
(painting procedure :1. Aluminium plate polished, 2. first 
coating , zinc chromate TT-P-1757 , 3. second coating , flat 
black paint , Part no. 8010 582 5382 ) in order to have 
emissivity close to unity, and was placed at a distance 
of 4m from the AGA (figure 6.1 shows the setup for the 
experiment). The water in the can was heated from 
ambient temperature up to its boiling point. The 
temperature of the can was measured with a 
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Figure 6.1 : Near field calibration setup 



thermocouple and every 2®C measurements of the thermal 
radiance were made using the AGA. 



This experiment 


was 


carried 


out 


twice and 


the 


"Isotherm Units" reading of 


the 


AGA for 


each 


temperature 


appears 


in 


Table 2 


. In 


figures 6.2 


and 


6.3 IU’s are plotted 


versus temperature 


for SW and 


LW 


respectively. 


In these 


same figures 


the calibration 


curves given 


by the 


manuf acturer 


were 


also plotted 


for 


comparison 


purposes 


. As 


can 


be seen the empirical 


calibration , 


data 


f alls 


pretty well into 


the 


manufacturers 


curves . 


The 


data 


of the curves noted 


"empirical 


calibration 2" 


is a 


little 


offset , it 


is 


believed , 


mainly 


because of 


inaccuracies in 


the 


thermocouple 


readings 


of 


the temperature ( offset 


of 


approximately 


3 to 4°C 


). 











C. FAR FIELD EXPERIMENT 



The thermal radiance of an aluminium plate of 1 

2 

m surface (lmXlm) was measured. The plate was painted 
black (Same painting procedure used as the one described for 
the Aluminium plate) in order to have emissivity close 
to unity. The experiment was carried out during 

two periods on the roof of Spanagel Hall (building 
232 , NPS ) . The plate heated by the sun can be assumed 
to have small fluctuations in its temperature throughout 
the experiment. Readings were taken with the plate at 
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TABLE 2 



NEAR FIELD MEASUREMENTS 



Temp. (°C) 


SW 


(ID) 


LW 


(IU) 


.st 

1 exp . 


2 n< * exp. 


i st 

1 exp . 


2 n< * exp. 


48 


51 


— 


77 


— 


50 


53 


- 


79 


- 


52 


57 


- 


81 


- 


54 


61 


- 


83 


- 


56 


63 


- 


85 


- 


58 


66 


- 


87 


- 


60 


71 


77 


90 


94 


62 


76 


83 


92 


97 


64 


80 


86 


93 


99 


66 


84 


90 


96 


101 


68 


89 


95 


98 


103 


70 


96 


99 


101 


105 


72 


99 


104 


103 


107 


74 


104 


109 


106 


109 


76 


106 


113 


107 


110 


78 


112 


122 


109 


114 


80 


116 


127 


111 


117 


82 


121 


134 


114 


120 


84 


130 


139 


116 


121 


86 


133 


148 


119 


124 


88 


142 


155 


121 


126 


90 


147 


165 


124 


132 


92 


155 


171 


127 


134 


94 


- 


177 


- 


136 


96 


- 


183 


- 


138 



Note: Ambient temperature : 23°C 
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SW CALIBRATION CURVES 
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Figure 6.2 : SW calibration curves 
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TEMPERATURE IN DEGREES CELCIUS 



TABLE 3 



FAR FIELD MEASUREMENTS 





SW 


(IU) 


LW 


(IU) 


Distance (ft) 


i si 

1 exp . 


2 nd exp. 


i st 

1 exp. 


0 nd 

2 exp . 


30 


48 


48 


76 


78 


60 


48 


48 


72 


78 


90 


46 


47 


68 


78 


120 


48 


44 


75 


75 


150 


47 


44 


74 


75 


180 


46 


44 


74 


75 


210 


44 


41 


74 


74 


240 


** 


41 


75 


73 


270 


** 


39 


74 


71 


300 


** 


37 


74 


71 


330 


** 


** 


75 


72 


360 


** 


** 


73 


74 


390 


** 


** 


71 


** 


Mean 


46.7 


43.3 


73.5 


74 


Standard 










deviation 


1.5 


3.7 


2.1 


2.3 


Corresponding 










temperature(®C ) 


43.5 


42 


47 


47 




* 


* 


* 


* 



Notes : * The actual temperature of the plate was 



in 


the 


l 3Z 


exp. 


: 42° C 


at 


the 


center , 










37°C 


at 


the 


corners . 


in 


the 


2 nd 


exp. 


: 44® C 


at 


the 


center 










38°C 


at 


the 


corners . 


Plate 


not 


detectable i 


by 


the 


systea. 
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different distance intervals (from the AGA) from 30ft 



up to 390ft away. 



The values 


obtained appear 


in 


Table 3 


. These 


values as it 


can be seen 


are lore or 


less the 


same , independent 


of distance. 








Specifically 


for the first 


far 


field 


set of 


measurements the 


plate temperature was 


42° C 


at the 


center and 37°. C 


at the corners 


and the 


average 


thermal values 


measured were 


46.7 


IU in 


SW (with 


standard deviation 


<r =1.5 IU) and 


73.5 


IU in 


LW (with 



or =2.1 10). These values correspond to 43.5 °C in SW 
and 47 °C in LW. 

For the second set of measurements the 



temperature 


of 


the 


plate 


was 


44° C at 


the 


center and 


38° 


C at 


the 


corners . 


The 


average 


thermal 


values 


measured were 


43.3 


IU in 


SW 


( or =3.7 


IU) 


and 


74 IU 


in 


LW ( <r 


= 2. 


3 IU) 


corresponding to 42° 


c 


and 


47° C 



respectively. 

D. CONCLUSIONS FROM EMPIRICAL CALIBRATION 

The conditions of the experiment ( Estimated 
accuracies : IU reading ± 5X , thermocouple reading ± 5X ) 

were not good enough to give good empirical calibration 
curves . 

The near field experiment showed that empirical 
calibration curves are not better than manufacturer’s 
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curves 



and the far field experiment showed that 



absorption of 


the atmosphere 


is Dot 


a 


big 


factor at 


small distances 


(up to 130m) 


• 








What has been shown is 


that 


for 


the 


purpose of 


the present 


work one 


should 


better 


use the 


manufacturer ' s 


calibration curves and 


keep 


in 


Bind that 


the results 


will have 


a tolerance 


of 


at least 


±5X ( IU reading 


accuracy ) . 
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VII 



SIGNATURE MEASUREMENTS 



A. MEASURED OBJECTS 

All the neasurenents, using the AGA , were made 

from the roof of Building 232 (Spanagel Hall) of the 

Naval Postgraduate School (NPS), (observers height 30m). 

The objects measured were in the vicinity of the 

NPS and specifically the following : 

(i) Different cloud formations at different ranges 
and altitudes. 

(ii) The following buildings: 

(a) The Monterey Sheraton Hotel located 

at 401 Alvarado St., Monterey CA. 93940. 

(b) Halligan Hall (building 234 of the NPS) 

(c) A warehouse located at 1260 9*"* 1 St., 

Monterey CA. 93940, viewed across the NPS. 

Figure 7.1 is a map of the Naval Postgraduate 

School showing the location of Halligan Hall and 

the warehouse at 9^"^ St. and figure 7.2 is a map of 

the city of Monterey showing the location of the 

Monterey Sheraton Hotel . 
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&IONTEREY BAY 




Figure 7.2 : Map of Monterey ,CA. 
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B. 



FORMULATION OF MEASUREMENTS 



An object will emit thermal radiation due to its 

temperature. The objective of the present work is to 
find this thermal radiance. 

The radiation received by the AGA is the sum of 
the object radiation, the radiation reflected by the 
object and the radiation emitted from the atmosphere 
( figure 7.3). 

If we make two assumptions : 

(i) That the observed object is opaque (ie. 

there is no transmission of IR radiation through the 
object ) , and, 

(ii) That the ambient temperature is uniform for 

all surroundings : 

we can assume that the three factors above are 

equivalent to three objects laid one on top of the 

other , with the same dimensions , that emit photon 

fluxes S , S and S , respectively, 
o a atm 

S q is the photon flux (photons/sec) , into the FOV 
of the receiver , emitted from a blackbody of 
temperature t Q at short distance, ie. no atmospheric 
absorption present . 

S is the photon flux (photons/sec), into the FOV 

a 

of the receiver , due to reflection of ambient 

radiation from the object . 



73 



ATMOSPHERE AGA 
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Figure 7.3 : The radiation budget. 



S is the photon flux (photons/sec) , into the 

FOV of the receiver , due to atmospheric emission and 
scattering. 

t is the object temperature. 

The first object , being the source itself , emits 
photon flux equal to S q , which attenuated via the 

atmosphere becomes 

object radiation = € r S 

o a o 

The second object , corresponding to the radiation 



reflected by 


the 


object 


, will emit photon flux equal 


to 


P S , 
o a 


where 


p is 

o 


the reflection 


coefficient of 


the 


source. 


From 


the 


assumption that 


the source is 



opaque we have that 




then the reflected photon flux is equal to 

(l-€ )S which attenuated through the atmosphere 

o a 

becomes 

reflected radiation = r (l-e )S 

a o a 

Finally the radiation emitted from the atmosphere 

will be, 

atmospheric radiation = ( 1 -t )S 

a atm 

Then we can write the radiation received by the 

AGA in the form of the following summation 
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(VII-1) 



S * =T € S + T (l-€ )S + (1-T ) S , 

o a o o a o a a at* 

=object radiation 

+ reflected radiation 

+ atmospheric radiation 

where , 



S f 
o 



a 



: is the radiation received via the 

atmosphere from the object surface 
(photons/sec) . 

: is the emissivity of the object . 

: is the atmospheric transmittance . 

The "AGA 780” has a linear photon counting 

detector. Thus the relationship between S values 

(photons/sec) and thermal values I (IU) is a linear 

one . 

I = C X S 

where C is a constant factor relating the signal 
in Isotherm Units to the received photon flux at the 
receiver (ie. C is an instrument response factor). 

This yields , 



I * = r € I + r (1~€ )I + (1-T )I 

o a o o a o a a atm 



( VII-2) 



where the I values are radiation detected by the 

AGA in IU and refer to the calibration curves of 

IU vs T °C provided by the manufacturer . 

Equation (VII-2) solved for I gives , 

o 
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( VII-3) 



I 

o 



€ T 

o a 





I 



atm 



I ’ is 

o 


the 


measured 


IU (thermal) 


value 


and I 

o 


is 


the thermal 


value 


of the 


object at 


the 


position 


of 


the object 


• By 


thermal 


value 


is 


meant 


the 


radiance of 


the 


object 


converted 


into 


IU through 



the calibration curves 

If we assune that the ambient temperature (t & ) 

is the same as the atmospheric temperature ( t . ) then 

atm 

the equivalent source thermal values from the 
surroundings and the atmospheric path will be equal, 
and equation (VII-3) can be rewritten as follows, 



I 

o 





1 ) I 



a 




I 

a 



( VII-4) 



(VII-5) 



Equation (VII-5) will be used to give the source 

thermal value of the object in IU , and from that 

2 

the radiance of the object (W/cm sr) will be computed. 

The procedure that was followed appears in the 
next section. 
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C. APPLICATION OF THE FORMULATION TO THE 
EXPERIMENTAL DATA 

(1) A Measurement with the AGA will provide a 

value I * (of irradiance equivalence) in IU. 

o 

(2) Entering the calibration curves of the AGA 
with the ambient temperature (t a ) we obtain the value 
I a (10) - 

(3) From the LOWTRAN 6 code with the meteorological 

and geographical data of each case we obtain the 

transmittance of the atmosphere . 

(4) The emissivity will be assumed to be 0.90 for 

clouds (in the SW) and 0.80 for buildings in the 

SW) . The LW emissivity is calculated assuming the 

emissivit ies in the SW are correct , as analyzed in 

the next chapter. 

(5) Using equation (VII-5) we obtain the thermal 

value Iq of the object . 

(6) From the thermal value I of the object we 

o 

find , using the calibration curves again , its 

temperature . The temperature of the object calculated 
from the SW data has been assumed to be the correct 

one ( see Chapter 8 ). 

(7) Using a TI-59 program (listing and analysis of 

that program appears in Appendix A) we compute the 

2 

in-band radiant flux density (W/cm ) of the object. 
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(8) If the observed object is a flat lambertian 



surface then the radiance , Lq 
be , 

€ M o 

L ql = (W/cm2 sr) 



from that object will 



( VII-6) 



The lambertian surface assumption was reasonable for 
most of the objects observed . 

For clouds where the thickness is much less 

than the range of observation another approximation 
might be considered ; the spherical surface approximation 
where the radiance of that object will be , 

€ M 0 

L = (W/cm2 sr) (VII-7) 

4n 



This 


last 


approximation 


was never 


used 


because 


in 


all the 


cases 


observed the 


extent of 


the 


clouds 


was 


greater 


than 


or about 


equal to 


the 


range 


of 



observation - 



D. SAMPLE CALCULATION 



(i) 


The 


radiance of 


a cloud in 


the SW 


and 


LW 


windows 


was 


measured 


The values 


obtained 


were 


20 



and 49 IU respectively. 

W = 20 IU 

I OLW ' = 49 IU 
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(2) Entering 



the 



calibration curves 



with the 



ambient 


temperature 


t 

a 


= 20°C 


we obtained the value 


I a’ 




^sw = 19 * 5 


IU 










I aLW = 47 ‘ 5 


IU 








(3) 


From the 10WTRAN 


6 code 


we obtained 






r = 0.1167 
a 


in 


the SW 


, and 






r = 0.3721 
a 


in 


the LW 






(4) 


Since we dealt with a cloud, its emissivity in 


the 


SW 


was assumed 


to 


be 0.90 


and the LW emissivity 


was 


calculated as 0. 


733 


(this 


calculation appears in 


the 


next chapter) . 












€ sw " 0,9 
€ lw = °* 733 










(5) 


Using equation 


(VII-5) we 


obtained 



I OSW 



20-19.5 

+ 19.5 

0.9X0.1167 



IU 



° r ’ I 0SW = 24 - 3 IU 



(6) From the 


calibration 


curves 


we found the 


temperature of 


the cloud 


as 


26°C 


using the SW I q 


value and 23°C 


using the LW 


I 

o 


value 


. This 


deviat ion 


is attributed 


to 


departure 


of 


the 


spectral 


radiance 


from the blackbody 


curve in 


the 


vicinity of 


the peak. 
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as found by E.E. Bell et al. [Ref. 14]. This "short 
wave" tenperature is thus assuied to be appropriate. 

(7) Using the TI-59 program the in-band radiant 
flux density was computed for the two wavebands. 

M ogw = 0.001178 W/cm 2 

Mqlw = 0.01699 W/cm 2 



(8) Finally 


dividing by 


n 


and 


multiplying 


by the 


corresponding 


emissivity 


we 


obtained the 


in-band 


radiance of the 


cloud in 


the 


SW 


and the LW 


windows . 



_!_”osw_ 

osw 

n 

or, 

0.9 X 0.001178 

L 0SW = = 3 • 375 x 10 ^ W/cn sr 



and, 



L 0LW 



!_”olw 

n 



or, 



0.733 X 0.0165 



'0LW 



= 3.964 



X 10 3 W/cm 2 sr 



Here the LW 
wave" emissivity 



radiance is computed using 
€ q as computed in Chapter 



the "long 
eight . 
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VIII. 



DATA ANALYSIS, PRESENTATION OF RESULTS 



A. DATA ANALYSIS 

1 . Emissivity 

The major uncertainty in the determination of 

cloud radiance from direct radiometric measurements is 
the determination of the appropriate value of 
emissivity . In the absence of a reference source at 

or near the location of the cloud this can only be 
inferred indirectly. 

The procedure adopted to determine source 

temperature was carried out using an assumed value 

(appropriate to literature values) of 0.9 for both 

long and short wavebands. In all cases the temperature 
of the object calculated from the SW data , assuming 

that the emissivity was 0.90 in both SW and 
LW (graybody approximation) , was greater than the 

temperature calculated using the LW data. 

An interpretation of this result is clear if we 

observe figures 3.2 and 3.3 . There we can see that 

for about the same conditions of measurements as the 
ones we had (cirrus-cumulus clouds, elevation angle * 15°) 
the cloud spectrum approaches the blackbody curves in 
the 3-5.6 pm range but is far different (less) in the 



82 



8-14 pm window . The same 
city scene in figure 3.16. 

This phenomenon is not 
transmittance because as it 
the transmittance in the SW 



was 


also 


observed 


for a 


due 


to 


the 


atmospheric 


can 


be 


seen 


in 


Table 5 


is 


less 


than 


that 


in the 



LW. 

The reasonable interpretation of this is that the 

emissivity in the SW is close to one (according to 

figures 3.2 and 3.3 ) , but that this is not the 

case in the LW . 

To account for this variation the emissivity in 

the SW has been assumed to be 0.9 for clouds 

and 0.8 for buildings . Using these values a 

temperature of the object is deduced . This temperature 
has been assumed to be correct , and from this 

tracking back a consistent emissivity for the the LW 

has been calculated. 

The procedure used is described below. 

2 . Calculation of the emissivity in the LW 

Equation (VII-5) solved for yields 



€ - 
O 



I ’ " I 
_o a_ 

r (I -I ) 
a o a 



( VIII-1) 



In this equation I * is the measured thermal 

o 

value (IU) of the object , I is the atmosphere 

a 

thermal value (IU) , I is the object thermal value 

o 
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taken 


from the 


calibration 


curves 


where we 


enter 


with 


the 


temperature of the 


object 


as computed 


with 


the 


SW 


data , 


and r is the 

a 


transmittance calculated 


with 


the 


LOWTRAN 


6 code. 









The following is a sample calculation of the 
emissivity in the LW band. 

In section VII. D the temperature of the cloud 
using the SW data was calculated to be 26°C. 

With this value entering the calibration curves 
we get, 

I = 53 IU 
o 

The other variables on the right hand side of 
equation VIII-1 are the same as in section VII. D, 

I *= 49 IU 
o 

I = 47.5 IU 
a 

r = 0.3721 
a 

Equation VIII-1 then gives an emissivity for the 

LW, 



to 

the 



€ oLW = °' 733 
This enissivi 

calculate the 



ty is 
in-band 



then used 
radiance 



LW. 



in 


step 


VII. D. (8) 


of 


the 


cloud in 



3 . Comments 

In some cases the computed emissivity in the LW 
did not fall in the permissible range ( 0<€<1 ). 
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This lay have been caused by a number of factors, 
of which the most likely are the following : 

(a) Invalid assumption that the emissivity is 0.9 

or 0.8 (clouds or buildings) in the SW. 

(b) Different transmittance from the one computed 

by the LOWTRAN 6 code . 

(c) Invalid assumption that the objects measured are 

opaque. 

(d) Error in the measurement. 

The analysis showed that the emissivity of the 

clouds in the SW is more or less close to unity; 

therefore the assumption of their emissivity being .9 
should be correct with an estimated error of ± .1 . 

The same holds for the emissivity of the buildings . 

The transmittance computed by the LOWTRAN 6 code 

might be in error , since the atmospheric model 

used (Midlatitude Winter ) might not match with the 

actual atmospheric conditions. No local meteorological 
measurements were made during the experiments. The 

magnitude of the error inserted in the computation of 

the emissivity , due to a transmittance error of 
a factor of 0.3 , will be a factor of three . 

The assumption of opacity of the objects is valid 
for the buildings and for clouds thicker than 60 m 
(fig. 3.11). 
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The Bain source of errors cones from the 
Measurement itself • An accuracy in the measured thermal 
value I q * of + 3 ID (which is about ±6%) could yield 

an enissivity error of about ± 0.9 which will cause 

the calculated enissivity to exceed the accepted range, 

( 0 <€< 1 ) . 

For those cases in which the computed enissivity 
in the LW was outside the valid range it has been 
assumed to be 0.7 for computation of radiance. 

B. PRESENTATION OF RESULTS 

The data collected appear in Table 4 along 

with information about the range of the object. Data 
on clouds was obtained from Monterey Air Traffic 

Control Tower ( MATCT ) and from Fritzsche Field 

( Fort Ord ) 

Table 5 shows the LOWTRAN atmospheric transmittance 

for each case along with the clear sky radiance. 

Table 6 shows the calculated radiance of the 

different objects along with their emissivity . 

Finally the average radiance values obtained appear 
in Table 7. 
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TABLE 4 



DATA COLLECTED 



Nuaber of 
Deasureient 


Object 


Meas « 

sw 


. (IU) 

LW 


Range 


i. 


SHERATON 


25 


52 


3.7km 


2. 


HALLIGAN 


26 


54 


100m 


3. 


WAREHOUSE 


27 


54 


150m 


4. 


CLOUD 


20 


49 


10km 


5. 


SHERATON 


20 


45 


3.7km 


6. 


CLOUD 


16 


44 


10km 


7. 


CLOUD 


15 


42 


10km 


8. 


CLOUD 


16 


43 


15km 


9. 


CLOUD 


15 


40 


15km 


10. 


HALLIGAN 


25 


50 


100m 


11. 


WAREHOUSE 


23 


50 


150m 


12. 


CLOUD 


16 


40 


10km 


13. 


CLOUD 


16 


43 


15km 


14. 


CLOUD 


17 


42 


15km 


15. 


CLOUD 


18 


44 


10km 


16. 


CLOUD 


15 


42 


10km 


17. 


LOW CLOUD 


21 


44 


7km 


18. 


LOW CLOUD 


21 


45 


7km 


19. 


CLOUD 


16 


44 


10km 



Note: (1) The ambient temperature for measurements 1 

through 4 was 20°C , and for measurements 5 
through 19 , 14°C. 

(2) The range to the clouds was estimated 
from weather information provided by the 
Monterey Air Traffic Control Tower ( MATCT) on 
the cloud coverage of the Monterey Bay area. 
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TABLE 5 



LOWTRAN OUTPUT 



Nuaber 

of 

aeas . 


Transmittance 


Radiance (W/cm sr) ! 


SW 


LW 


SW (X 10 4 ) 


LW (X 10 3 ) 


i. 


0.4296 


0.7574 


0.9966 


0.8437 


2. 


0.7954 


0.9618 


0.4997 


0.1593 


3. 


0.7954 


0.9618 


0.4997 


0.1593 


4. 


0.1167 


0.3721 


0.944 


1.866 


5. 


0.4296 


0.7574 


0.9966 


0.8437 


6. 


0.1067 


0.4972 


1.231 


1.547 


7. 


0.2065 


0.5157 


1.021 


1.497 


8. 


0.1247 


0.3872 


0.9753 


1.832 


9. 


0.1247 


0.3872 


0.9753 


1.832 


10. 


0.7954 


0.9618 


0.4997 


0.1593 


11. 


0.7954 


0.9618 


0.4997 


0. 1593 


12. 


0.1067 


0.4972 


1.231 


1.547 


13. 


0.1247 


0.3872 


0.9753 


1.832 


14. 


0.1247 


0.3872 


0.9753 


1.832 


15. 


0. 1935 


0.4972 


1.001 


1.547 


16. 


0.1935 


0.4972 


1.001 


1.547 


17. 


0.2243 


0.5403 


1.014 


1.45 


18. 


0.2243 


0.5403 


1.014 


1.45 


19. 


0.1935 


0.4972 


1.001 


1.547 



The above values are path average transmittance 
and path radiance computed using the LOWTRAN 6 
code. 
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TABLE 6 



COMPUTED RADIANCE 



Nunber 

of 

meas . 


SW eaiss. 


SW Radiance 
W/ ci sc 
X10 H 


LW eaiss. 


LW Radiance 
W/cin sr 
X10 J 


1 . 


0.8 


3.888 


0.838 


5.108 


2. 


0.8 


3.5345 


0.778 


4.539 


3. 


0.8 


3.649 


0.711 


4.210 


4. 


0.9 


3.375 


0.733 


3.964 


5. 


0.8 


3.476 


0.985 


5.703 


6. 


0.9 


2.951 


0.7 ** 


3.558 


7. 


0.9 


1.928 


0.7 *** 


2.923 


8. 


0.9 


2.851 


0.32 


1.601 


9. 


0.9 


1.791 


0.7 *** 


2.825 


10. 


0.8 


3.769 


0.601 


3.612 


11. 


0.8 


3.206 


0.742 


4.138 


12. 


0.9 


2.951 


0.7 ** 


3.558 


13. 


0.9 


2.851 


1 


5.004 


14. 


0.9 


3.879 


0.49 


2.823 


15. 


0.9 


4.105 


0.7 ** 


4. 144 


16. 


0.9 


1.859 


1 


4.107 


17. 


0.9 


5.767 


0.22 


1.525 


18. 


0.9 


5.767 


0.516 


3.576 


19. 


0.9 


2.851 


1 


5.004 



The emissivity in the SW band is assumed as 0.9 

for clouds and 0.8 for buildings . 

The emissivit ies marked ** or *** are assumed 0.7 

because the calculation yielded emissivity out of 
accepted limits ( >1 and <0 respectively ) 
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TABLE 7 



AVERAGE RADIANCE VALUES 



1 

1 


SW RADIANCE 1 


LW RADIANCE 


1 

1 

1 


2 -4 ' 

W/cn sr X 10 1 


W/ci^sr X 10 ^ 


1 

1 

CLOUD : 

FORMATIONS ! 

I 

l 


i 

i 

3.302 : 

1 

1 

1 

1 


3.544 


l 

l 

MONTEREY ! 

SHERATON ! 

I 

l 


1 

1 

3.682 ! 

i 

• 

i 

i 


5.406 


i 

l 

HALLIGAN \ 

HALL ! 

1 

i 


i 

i 

3.652 I 

• 

i 

i 

i 


4.076 


i 

i 

WAREHOUSE ! 
AT 9th ST. 1 

i 

i 


i 

i 

3.428 ! 

i 

i 

l 

i 


4.174 
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IX . CONCLUSIONS / RECOMMENDATIONS 



A. CONCLUSIONS 

The radiance leasureaents obtained appear 

analytically in chapter 8, 

A comparison of radiance from a cloud, for 
example, and the radiance of the background sky, shows 

that the cloud would always have, in these cases , a 

positive contrast with respect to its background. This 
shows the necessity of the existence of more than a 

background supression scheme ie. a system that, after 

the subtraction of the background radiance, will 

differentiate between real targets and clutter. 

B. RECOMMENDATIONS 

During the present work the main problem 

encountered was the translation of IU into meaningful 
radiance values. The problem was solved by using very 

general approximations ; the Lambertian surface 

approximation and the spherical surface approximation. 



For better 


results a 


cal ibrat ion 


curve 


of IU 


versus radiance 


units 


would be 


useful 


The 


representative 


of 


the AGA 


Corporation 


(Mr. 


Jack M. 


Patterson , 108 


Arena 


st . , El 


Segundo , CA. 


90245 


, tel . 


(213)- 3226257) 


has 


been 


contacted on 


that 


subject 
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and a telex has been sent to the nanuf acturing 

coipany in Sweden. 

At the time of completion of this work no answer 
yet has been received on the matter. 

It is recommended that this subject be pursued 



because , 


for any application, 


and 


specifically 


for 


radiance 


measurements 


like 


the 


ones 


made here , 


it 


will be 


very helpful 


to 


have 


on 


hand a calibration 



2 

of IU measured versus radiance (W/m sr) emitted by an 
object . 

An interesting experiment will be to measure 

simultaneously an object with both the AGA and the 
IRSTD and obtain IRSTO signal calibration versus 

radiance . 
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APPENDIX A 



TI-59 PROGRAM FOR CALCULATION OF THE IN-BAND FLUX 

A. GENERAL 

The TI-59 program used to calculate in-band flux 

was originally written by Robert Pitlak (Appollo Lasers, 
6357 Arizona Circle , Los Angeles CA 90045 ) for an HP-67 
calculator. 

The inputs of the program are: 

(1) The initial and final wavelengths that define 
the band . 

(2) The temperature of the object (in °K, or in 
®C , or in ° F) . 

(3) The number of intervals in the numerical 

integration ( this no. must be even ). 

2 

The result may be given either in Photons/cm sec 

2 

or in W/cm 

B. THEORY 

Plank’s blackbody radiation law may be written in 
the form: 





or 




(q/X) dX 
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where 



and 



X 4 (e x -l) 



c = 1.883651556 X 10 
o 



X = 



_ 2 _ 

XT 



23 



c L = 37418.42875 
c 2 = 14388.32334 





The 


above integrations are 


carried out 


by the 


TI- 


59 


numerical ly 


in 


any (even) no. of 


intervals 


specif ied 


by the 


user • 










c. 


PROCEDURE 














To use the program 












(i) 


Load the 


card 


into the TI-59 (both faces 


to 






be read) 


• 












(2) 


Press RST . 












(3) 


Enter X^ 


in i 


microns , 


then press 


R/S. 






(4) 


Enter X 2 


in i 


microns , 


then press 


R/S. 






(5) 


Enter T 


( temperature) 


then press 










A 


if in 


° K 












B 


if in 


°F 












C 


if in 


°C 










(6) 


Enter 


n : 


number 


of intervals 


of 


the 










numerical 


integration . 














( n MUST 


be EVEN ) 








(7) 


Press 


D for 


answer 


2 

in W/cm 










or 


E for 


answer 


2 

in Photons/sec cm 
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D. LISTING OF THE PROGRAM 



listing 


of 


the 






* 


055 








056 








057 








053 


000 


42 


STD 


059 


001 


01 


01 


060 


002 


99 


PRT 


061 


003 


91 


R/S 


062 


004 


42 


STD 


063 


005 


02 


02 


064 


006 


99 


PRT i 


065 


007 


91 


R/S 


066 


003 


7 6 


LBL 


067 


009 


11 


ft 


0 6 3 


Old 


42 


STD 


069 


01 1 


00 


00 4 


070 


012 


02 




071 


013 


06 


3 


072 


014 


00 


0 


073 


015 


00 


0 


074 


016, 


66 


0 


075 


017 


00 


0 


0 ■ 6 


013 


69 


DP 


077 


019 


04 


04 


073 


020 


43 


RCL 


079 


021 


00 


00 


030 


022 


69 


DP 


031 


023 


06 


06 


032 


024 


69 


DP 


033 


025 


00 


00 


034 


026 


91 


R/S 


035 


027 


7 6 


LBL 


036 


023 


12 


B 


037 


029 


42 


STD 


033 


030 


00 


00 


039 


031 


02 




090 


032 


01 


T 


091 


033 


00 


o 


092 


034 


00 


0 I 


093 


035 


00 


0 


094 


036 


00 


0 


095 


037 


69 


DP 


096 


033 


U4 


IJ4 i 


097 


039 


43 


RCL 


093 


040 


00 


00 


099 


041 


69 


DP 


100 


042 


06 


06 


101 


043 


69 


DP 


102 


044 


00 


00 


103 


045 


75 


- 


104 


046 


03 


3 


105 


047 


02 


2 


106 


043 


95 


= 


! 107 


049 


65 


V 


103 


0 5 0 


05 


5 


109 


051 


55 




110 


052 


09 


y 


111 


053 


35 


+ 


112 


054 


02 


2 


, 113 



follows : 








114 


02 


2 


174 




LMX 


115 


03 


•-* 


175 


42 


STD 


116 


03 




176 


06 


06 


117 


04 


4 


177 


36 


PGM 


113 


55 


- 


173 


09 


09 


119 


43 


RCL 


179 


14 


D 


120 


07 


07 


130 


22 


I MV 


121 


55 


-r 


131 


52 


EE 


122 


43 


RCL 


132 


99 


PRT 


*. 123 


00 


00 


133 


93 


ADV 


I 124 


54 


> 


134 


91 


R/S 


J 125 


54 


) 


135 


7 6 


LBL 


; 126 


75 


- 


136 


15 


E 


127 


01 


1 


137 


93 


ADV 


123 


54 


> 


133 


42 


STD 


129 


54 


i 


139 


05 


05 


130 


35 


1/X 


190 


99 


PRT 


131 


65 


X 


191 


93 


ADV 


132 


43 


RCL 


192 


53 


< 


133 


03 


03 


193 


43 


RCL 


134 


54 


) 


194 


02 


02 


135 


92 


RTN 


195 


■ 5 


- 


136 


7 6 


LBL 


196 


43 


RCL 


137 


14 


D 


197 


01 


01 


133 


93 


ADV 


198 


54 


> 


139 


42 


STD 


199 


55 


• 


140 


05 


05 


200 


43 


RCL 


5 141 


99 


PRT 


201 


05 


05 


i 142 


93 


ADV 


202 


95 


= 


1 143 


53 


< 


203 


42 


STD 


144 


43 


RCL 


204 


03 


03 


145 


02 


02 


* 205 


01 


1 


146 


75 


- 


. 206 


93 




: 147 


43 


RCL 


207 


03 


3 


143 


01 


01 


203 


03 


3 


149 


54 


) 


209 


03 


3 


150 


55 


- 


210 


06 


6 


151 


43 


RCL 


211 


05 


5 


152 


05 


05 


s 212 


01 


1 


153 


95 


= 


, 213 


06 


6 


! 154 


42 


STD 


i 214 


52 


EE 


155 


03 


03 


215 


02 




156 


03 


•-s 

O 


• 216 


03 


3 


157 


93 


• 


. 217 


42 


STD 


153 


07 


\ 


w 1 'J 


03 


03 


159 


04 


4 


' 219 


04 


4 


160 


01 


1 


220 


42 


STD 


161 


05 


5 


. 221 


09 


09 


162 


00 


0 


i 222 


01 


1 


163 


07 


7 


1 223 


•to 


I MV 


164 


04 


4 


224 




LMX 


165 




EE 


! 225 


42 


STD 


1 66 


04 


4 


226 


06 


06 


1 67 


42 


STQ | 


227 


36 


PGM 


163 


03 


03 j 


223 


09 


09 


169 


05 


5 


229 


14 


D 


170 


42 


STD 


230 


•to 


IMV 


171 


09 


09 


231 


52 


EE 


172 


01 


i ! 


232 


99 


PRT 


173 


22 


I MV ! 


233 


93 


ADV 








234 


91 


R/S 



program 

07 7 

03 3 

35 = 

42 STD 

00 00 

91 R/S 

76 LBL 

13 C 

42 STD 

00 00 

oi i 

05 5 

00 0 

00 0 

00 0 

00 o 

69 OP 

04 04 

43 RCL 

00 00 

69 DP 

06 06 

69 DP 

00 00 

43 RCL 

00 00 

35 + 

02 2 

07 7 

03 3 

95 = 

42 STD 

00 00 

91 R/S 

76 LBL 

16 A' 

42 STD 

07 07 

53 < 

53 < 

43 RCL 

07 07 

45 Y* 

43 RCL 

09 09 

65 x 

53 < 

53 < 

43 RCL 

06 06 

45 YX 

53 < 

01 1 

04 4 

03 3 

03 3 

03 3 

93 . 

03 3 
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